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^ : ABSTRACT 
O , 

Many gamma-ray burst counterpart searches are being conducted in the 
^ _ optical-UV band. To both predict detectabihty and understand the meaning 

Tjj- ' of any detections or upper hmits, we extrapolate gamma-ray spectra from 54 

in ■ 
o 

■ concurrent with gamma-ray emission and do not consider quiescent or fading 



O 



counterparts. We find that the spectrum must be steeper (greater flux at low 



^ ■ energy) than a simple extrapolation of the gamma-ray spectrum for more than 

1^ • one simultaneous optical flash to be observable per year by current searches. 

o 

c/3 . Subject headings: gamma rays: bursts 

> 

tj ■ 1. Introduction 

P.. 

Despite years of study and volumes of data, gamma-ray bursts (GRBs) remain among 
the most mysterious cosmic phenomena. Perhaps the greatest impediment to theoretical 
progress is that the source of GRBs is unknown. Although it is widely thought that neutron 
stars are involved, there is no consensus on even the broadest details of the physical processes 
involved. The main reason for this is that the distance to GRBs is not known. Many have 
interpreted the observations of GRB isotropy and spatial inhomogeneity (Fishman et al. 
1994) by the Compton Gamma Ray Observatory'' s Burst and Transient Source Experiment 
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(BATSE) as indicating that bursts are at cosmological distances. Attempts have been 
made to verify this conclusion with BATSE data (Davis et al. 1994; Norris et al. 1994) but 
the claims are controversial (Band 1994; Brainerd 1994) and the possibility of a Galactic 
origin has not been ruled out (c.f., Podsiadlowski, Rees, & Ruderman 1995). Uncertainty 
in the distance to GRBs leads to uncertainty in the energy required to generate observed 
GRBs and hence, confusion about the source physics. If the distance scale to GRBs can 
be determined, it will be the greatest step yet taken toward understanding these enigmatic 
events. 

An excellent way to establish the GRB distance scale is to discover counterparts at 
other wavelengths. Counterparts could also reveal valuable information about burst sources 
and environments which would have an enormous impact on our physical understanding 
of these events. For example, if GRBs were found only in a particular environment, then 
properties unique to that environment are likely to be responsible for the burst. Many 
groups (see reviews by Schaefer 1994 and Greiner 1995) have attempted to find counterparts 
which flare along with the burst, are observable in quiescence, or appear as afterglows. To 
date, convincing counterparts have been found only for soft gamma repeaters, which many 
believe to be unrelated to GRBs (although this opinion is not universal; see Rothschild & 
Lingenfelter 1996). Several theories have been been proposed which predict emission at 
wavelengths other than gamma-rays (a number of them are listed in Schaefer 1994) although 
the guidance they provide is slim considering that many were advanced at a time when 
bursts were thought to originate at distances ^ 1 kpc, much closer than is now believed. 

Since we have no clear expectations as to the appearance of a burst source in quiescence 
or immediately after the gamma-ray event, we assume in this work that all optical emission 
occurs at the same time as the high energy emission and is the low energy tail of the 
gamma ray spectrum. Just as we extrapolate observed X-ray and gamma-ray emission to 
the optical-UV band, some model high energy spectra can be extrapolated to lower energy. 
Katz (1994) and Tavani (1996a,b) proposed a model where particles in a relativistic shocked 
plasma with an equilibrium energy distribution radiate a synchrotron spectrum which can 
extend to optical energies. Meszaros & Rees (1993) modeled the multiwavelength emissions 
of a fireball which interacts with the surrounding medium, resulting in a reverse shock, 
with emission from a number of regions. A variety of different high energy and optical-UV 
spectra are possible, depending on the physical parameters (e.g., densities and magnetic 
field) ; for some parameters an optical synchrotron spectrum similar to the Katz and Tavani 
models is emitted, while for others the optical emission can be greater, or self-absorbed. 

Several models predict concurrent low energy emission based on the reprocessing of 
gamma-rays in a neutron star magnetosphere (Hartmann, Woosley, Sz Arons 1988), an 
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accretion disk (Epstein 1985) or a stellar companion (London & Cominsky 1983). These 
models all assume nearby (~ 100 pc) somxcs and consequently much smaller source 
energies 10^^ ergs) than is required for bursts in the Galactic halo or at cosmological 
distances lO'*^ ergs); it is unclear what effect such a large energy release will have on the 
reprocessing region. While these reprocessing models predict optical re-emission while the 
burst is still in progress, the re-emitted flux should be proportional to the total gamma ray 
energy flux and not to an extrapolation of the gamma ray spectrum. 

Lacking compelling theoretical expectations, we take an empirical approach to 
determine how bright simultaneous counterpart emission may be by extrapolating spectra 
observed by BATSE to optical and ultraviolet energies (§2). Several experiments capable of 
detecting simultaneous optical emission from a GRB are described (§3) and the implications 
of our results for these experiments are discussed (§4). 

2. Analysis 

We base our extrapolated fluxes on spectral fits to 54 bright GRBs from the first 13 
months of BATSE operation (Band et al. 1993). The bursts were selected with a brightness 
criterion using the peak count rate in the 50 — 300 keV energy band. Most bursts which 
met the intensity requirement were included in the sample, but there were a few exceptions. 
Bright bursts shorter than Is were left out of the sample while some dim bursts with 
interesting time structure were included; the number in both cases was small. Although the 
sample is not statistically complete, it is representative of bright bursts. 

In this work, we ignore interstellar absorption. While this effect will render sources 
which must be observed through the Galactic plane invisible in the optical-UV band, it 

is far less important above the plane. Near the Galactic poles, absorption typically dims 
sources by Amy ~ 0.3 (Mihalas & Binney 1981, p. 181). The exact amount of absorption 
depends on the amount of interstellar matter along the line of sight, something which is 
not a simple function of elevation angle above the plane because of the clumpy nature 
of the interstellar medium. Given the large range of magnitudes which results from our 
extrapolation of the spectrum over four energy decades under very simple assumptions, the 
additional complexity of modeling interstellar absorption of Amy ~ 0.3 — 1 off the Galactic 
plane is not warranted. 

Band ct al. (1993) fitted spectra averaged over the entire burst to a four parameter 
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function of the form 



VkeV-s-cm / 



where A, a, f3, and Eq were fitted to the observed spectra and A' was chosen to make 
the function continuously differentiable everywhere. The energy range covered by the 
fit varied, with typical lower and upper energy cutoffs of 10 — 30 keV and 1 — 3 MeV, 
respectively. This is a differential photon number spectrum which must be multiplied by 
energy, Fe{E) = Ne{E)E, to get the energy flux. Note that the spectral indices in eq. (1) 
do not have an explicit minus sign. 

Here we extrapolate the fitted spectra down to 1 keV where they break to a power law 
with a different spectral index which describes the optical through soft X-ray spectra. The 
index of the low energy extrapolation is varied from —4 to 4 in integer steps. In addition, 
we include an energy index of 1/3, which is physically interesting, as we discuss below. 
Finally, a low energy power law spectrum of index a + 1, an unmodified continuation of the 
gamma ray spectrum, is extrapolated to the optical range. 

Discussions of counterpart emission are conventionally based on burst fluence; expected 
magnitudes are presented as if all the emission occurred in one second (Schaefer 1981). 
Although this is not consistent with the assumption of strictly simultaneous optical 
emission, it is a useful exposure-time independent standardization. We calculate fiuence 
based magnitudes in the energy bands of the U, B, and V filters and the 5 — 7 eV bandpass 
of the UV camera array on the High Energy Transient Explorer (HETE, discussed in §3). 
For the first three bands, we convolved spectra with the filter response functions (Budding 
1993, p. 46) and converted the fiuxes to magnitudes. The i^ETi? bandpass is not a standard 
wavelength band so we integrate the flux over the 5 — 7 eV band. Expected magnitudes 
for a one second exposure are given in Table 1 and distributions for the a + 1 and 1/3 
extrapolations are shown in Figures 1 and 2 (the distributions for the integral indices 
are identical to the distribution for 1/3, but with the magnitude axis shifted). The true 
differential distribution is a convolution of the intensity distribution and the distribution of 
spectral shapes. Since the burst intensity distribution increases with decreasing intensity 
(perhaps until the spatial inhomogeneity becomes particularly severe), we expect the 
optical magnitude distribution to increase with magnitude. However, our calculated optical 
distributions are based on a sample of bright bursts, and therefore the distributions are 
truncated on the large magnitude (low optical flux) side. If the optical searches can 
reach magnitudes larger than or comparable to the peaks of the distributions in Figures 1 
and 2, then the number of bursts which these searches may detect will be underestimated. 
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However, as will be demonstrated below, this is not an issue for physically interesting 
extrapolations. 

Although fluence-based measures of brightness are the standard, flux-based measures 
are more meaningful when detector exposure durations are comparable to the burst 
duration (note that we assume the optical emission has the same duration as the gamma-ray 
emission). This distinction is unimportant when the exposure time is longer than the burst 
since the entire burst fluence will be observed. Indeed, the search for optical transients 
began with the inspection of photographic plates which had exposure times much longer 
than the typical burst duration. If, however, the exposure time is shorter than the 
burst, then only a fraction of the total flux is observed and fluence-based measures will 
overestimate the brightness of a burst. 

Another concern when dealing with real experiments is spectral response. All 
the instruments we consider except HETE are ground based and sensitive between 
~ 4000 — 7000A. Although this energy band is very wide, the convention is to define 
sensitivities in terms of V magnitudes. Therefore, we convolve spectra with the response 
function of the V-filter and convert the flux to magnitude for comparison. This will 
introduce slight errors in the comparison of observability because the GRB spectra are 
assumed to be power laws while the sensitivities of the instruments were determined using 
stars, which have black body continua. These errors are not large however and a proper 
treatment would require detailed knowledge of the complicated response functions of the 
various experiments. 

For each instrument we calculated the expected observed flux within each experiment's 
exposure time. In all cases, the detector response was assumed to be linear (no saturation). 
If a burst was shorter than the exposure time, the entire flux of the burst was averaged over 
the exposure time (equivalent to the fluence-based magnitude with dilution for exposure 
times > Is). If the burst was longer, we used the average flux over the exposure time. This 
process can be expressed mathematically as 



where fv is the V band extrapolated flux and /o is the V band flux of a magnitude zero 
object (3.08x10-*^ ergs/cm^-s— Budding 1993, p. 46). T is defined as 



where t^xp is the instrument exposure time and tdur is the burst duration. The case 
texp < tdur IS Only an approximation to the true instantaneous fiux because intensity and 



my = 2.51og(r/y//o). 



(2) 




(3) 
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Fig. 1. — The number of bursts at a given fluence-based magnitude (entire burst occurs in 
one second) assuming the flux extrapolates as a power law of index a + 1 from the gamma- 
rays and a one second exposure. Panels (a)-(d) are the V, B, U, and 5 — 7 eV [HETE) 
bands, respectively. Note: Only bright bursts are included. If all BATSE bursts had been 
included in the sample, this differential distribution would continue to rise as the magnitude 
increased, not peak as shown here. 
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Fig. 2. — Same as Figure 1, but for an energy index of 1/3. 
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spectra vary within a burst (Ford et al. 1995). Although the spectral changes are not 
as radical as intensity fluctuations, spectral variations are amplified by the four decade 
extrapolation from the lowest energy gamma-rays included in the original spectral fits. 

Better sensitivity can be obtained for the short exposure instruments by adding 
subsequent images together. We did not consider this effect since approximately one 
magnitude is gained for a typical burst. This boost in sensitivity is likely to be more than 
compensated by our neglect of slew times for telescopes which respond to a burst trigger 
since typically the brightest emission occurs early in the burst. 

3. Instrumentation 

Instruments which can detect simultaneous optical emission from GRBs can be grouped 
into two categories. The first has large fields-of-view and long exposure times (> 30 min). 
These telescopes were designed for meteor patrols, wide field imaging or variable star 
patrols, but their large ficlds-of-view make them well suited for GRB detection. Many 
of these instruments arc outlined in Greincr ct al. (1994). Those most likely to image a 
field containing a GRB are sky patrols from Sonncbcrg and Ondfejov (the Odessa and 
Dushanbe sky patrols have been suspended because of economic difficulties in Ukraine 
and Tadshikistan — J. Greiner 1996, private communication). Although the patrol plates 
have a variety of sensitivities and exposure times, we concentrate on sky patrol plates from 
Ondfejov because its field-of-view is particularly wide and it has three coincidences with 
bursts in the Band et al. (1993) sample (Greiner et al. 1994). 

The second class of instruments we consider is characterized by short exposure times. 
These instruments are much more sensitive to GRBs, which seldom last more than a 
few minutes (the longest burst in our sample is 130s). We consider four experiments 
here: the Explosive Transient Camera (ETC, Vanderspek et al. 1994), the Gamma- Ray 
Optical Counterpart Search Experiment (GROCSE, Akerlof et al. 1994), the Gamma-ray 
To Optical Transient Experiment (GTOTE, S. Barthelmy 1995, private communication), 
and the Ultraviolet Transient Camera Array on the HETE spacecraft (Ricker et al. 1992; 
Vanderspek et al. 1995). All of these instruments monitor a portion of the sky but the 
ETC, GROCSE, and GTOTE can slew to preliminary BATSE positions distributed on the 
BACODINE network (Barthelmy et al. 1994). BACODINE requires up to five seconds to 
localize a GRB and another few seconds to distribute the coordinates on the network. After 
that, the telescope must slew to the proper position. These delays will often result in bright 
portions of the burst being missed by the telescope. We ignore these concerns in this paper 
so our extrapolated brightness may be somewhat optimistic. Details of all the instruments 
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mentioned here are presented in Table 2. 

4. Discussion 

The flux-based magnitudes of bursts expected under our assumptions are given in 
Table 3. Figures 3a and b show the V-band brightness distribution of our sample for 
extrapolations with energy indices of o; + 1 (the low energy gamma-ray power law) and 
1/3, respectively, for the ETC. Although the ETC, GROCSE, and GTOTE have varying 
integration times, only six of the 54 bursts in our sample were shorter than five seconds (the 
exposure time of the ETC). These bursts are slightly brighter in GROCSE and GTOTE. 
The difference is only a few tenths of a magnitude though and the resulting distributions 
have the same maximum and median brightness. The overall distributions are almost 
identical as well so that Figure 3 applies to all three of these experiments. Figures 4a and b 
show similar distributions for HETE. 

An extrapolation with an energy power law index of 2 could correspond to the 
Rayleigh- Jeans portion of a thermal spectrum. For the Rayleigh-Jeans power law to extend 
to 1 keV the radiating plasma should have a temperature above ~ 1 keV, which is possible 
in the extreme conditions expected in the burst environment. The thermonuclear model of 
Woosley & Wallace (1982) predicted an X-ray tail with a temperature greater than 2 keV; 
this model applied to local bursts, and the predicted emission was not concurrent, but the 
model does show that high temperatures are plausible. However, Table 3 shows that if this 
is the only component of concurrent GRB emission at lower energies, GRBs will not be 
observable at optical or UV wavelengths (assuming no afterglow). 

A power law index of 1/3 might result from low energy synchrotron emission (Katz 
1994; Tavani 1996a,b). The low energy limit of the synchrotron spectrum emitted by 
a single electron is a power law with an index of 1/3 below a characteristic energy of 
Eg = {3/4'iT){ehB/'meC)'~f^, where B is the magnetic field, and 7 is the Lorentz factor 
corresponding to the electron's energy. The synchrotron power law indices of ~ — 1 
familiar from radio astronomy result from a power law electron distribution. If the electron 
distribution has a low energy cutoff 7c, then the low energy single electron synchrotron 
power law of index 1/3 will be observed below -E's(7c). Note that the energy spectral index 
(l+a) of the low energy component of the spectral function in eq. (1) is typically of order 
1/3. In Figures 2, 3b, and 4b we present the optical- UV burst distributions for an energy 
index of 1/3. For this case the brightest burst in our sample would have appeared as a 
my ~ 13 transient, just at the sensitivity limit of GROCSE. Therefore, we doubt current 
instruments can effectively test this model of counterpart emission. 
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Fig. 3. — The number of bursts at a given flux-based V magnitude for a 5 s exposure 
assuming tlie flux extrapolates as a power law of index a + 1 (panel a) or 1/3 (panel b) 
from the gamma-rays. The dashed hnes indicate the sensitivities of the ETC (my = 11) and 
GTOTE and GROCSE (my = 13) experiments. 
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Fig. 4. — The distributions of extrapolated fluxes for a 4s exposure in the HETE band of 
5-7 eV for energy power law index a + l (panel a) or 1/3 (panel b). The dashed line indicates 
HETEs sensitivity. 
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The results of previous optical searches rule out some of the spectral indices considered 
in this simple model. In particular, Greiner et al. (1994) examined sky patrol plates from 
several sites which were coincident with GRBs. Three plates from the Ondfejov patrol were 
coincident with bursts in our sample. These plates had four hour exposures with sensitivity 
to events brighter than my ~3. The three bursts which overlap with our sample (3B911127, 
3B920525, and 3B920530) would have produced optical images at least this bright if the 
spectrum was a power law of index less than —3 below 1 keV. (Pubhshed before the 3B 
catalog was released, this result is still valid since the difference in the positions of these 
bursts in the 2B and 3B catalogs is less than the plate dimensions — J. Greiner 1996, private 
communication) This rules out more steeply falling power laws, which is not a severe 
constraint. 

Lacking compelling physical guidance for our expectations, the simplest hypothesis is 
that spectra in the optical band have the same index as at higher energies (the a + 1 case) . 
If this is true, then the ETC, GROCSE, and GTOTE have a chance of detecting the very 
brightest bursts which trigger BATSE and arc distributed by the BACODINE network. 
The ETC would have seen only the brightest burst in the sample while GROCSE and 
GTOTE would have seen four. HETE would not fare as well, being unable to detect any 
GRBs in the UV if the a + 1 extrapolations hold. 

Being able to detect a GRB does not imply it will be observed. Assuming that each 
site can image half the sky half the time (night only), then there is a 25% chance a BATSE 
GRB can be slewed to. For simplicity we neglect other factors which reduce the duty cycle, 
such as inclement weather and moonlight. Therefore, GROCSE and GTOTE might be able 
to detect one burst per year for the a + 1 extrapolation and the ETC might record a positive 
detection every four years. All of these detectors have a monitor mode which could catch a 
burst BATSE does not see. BATSE has ~33% sky coverage, so the number of potentially 
detectable GRBs is tripled assuming a site can view all the sky all the time. However, 
75% of these are lost to a given ground based site when earth blockage and daytime are 
included. Most of the rest are inaccessible because the fields-of-view of these instruments is 
small. The ETC has the largest field-of-view and it covers only 10% of the sky accessible to 
it. Therefore, detection in the monitor mode of any of these instruments is unlikely for this 
extrapolation. 

The instruments discussed here would have difficulty observing simultaneous optical 
emission from a GRB if optical photons are generated as part of the same spectrum as the 
gamma-rays. This implies that if concurrent emission is observed more than about once 
a year, the energy flux must follow a power law of index greater than a; + 1, indicative of 
another spectral component since GRB spectra are concave down above 10 keV. However, 
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there are indications of emission in excess of simple extrapolations below ~ 10 keV in 10% 
of GRBs (Preece et al. 1996). Therefore, aside from the discovery of a GRB counterpart, 
the observation of simultaneous optical emission would indicate a multifaceted emission 
process, providing new and important information about the burst environment. 

5. Summary 

Using spectra of 54 bright GRBs from BATSE, we have determined the expected 
optical brightnesses for simultaneous optical emission for an empirical extrapolation of the 
observed gamma-ray spectrum. We find that in order for concurrent optical emission to 
be observable, the spectral index of a power law extrapolation must be negative at long 
wavelengths (more flux at lower energies). A simple extrapolation to optical wavelengths is 
only marginally detectable, providing about one observable burst per year for the GROCSE 
and GTOTE telescopes. This imphes that if optical transients concurrent with burst 
emission are observed more frequently, an additional spectral component must exist. 

We are indebted to S. Barthelmy and R. Vanderspek for sharing details of their 
instruments. We also thank R. Narayan for suggesting this project and W. Coburn, 
G. Huszar, D. Marsden, and L. Peterson for their assistance. We appreciate the insightful 
comments of the referee, J. Greiner. This work was supported by NASA contract 
NAS8-36081. 
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Table 1. Fluence Based Expected Brightness^ 
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max(,/ViM-)'' 


iii(HLiaii( /7 iLi:)^' 


a + 1 
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15.3 


8.0 


15.4 


7.2 


14.4 


1.03x10"" 


1.72x10-" 


-4 


-17.3 


-12.3 
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-11.1 


-16.0 


-11.1 


3.13x10-3 


3.57x10-5 


-3 


-10.7 


-5.7 


-9.7 


-4.7 


-10.0 


-5.0 


1.81x10-5 


2.06x10-^" 


-2 


-4.1 


0.9 


-3.4 


1.6 


-3.8 


1.2 


l.OSxlO-'' 


1.20x10-8 


-1 


2.5 


7.5 


3.0 


8.0 


2.3 


7.3 


6.22x10-1" 


7.08x10-12 





9.1 


14.1 


9.3 


14.3 


8.4 


13.4 


3.70x10^12 


4.21x10-" 


1/3 


11.3 


16.3 


11.4 


16.4 


10.4 


15.4 


6.71x10-13 


7.64x10-15 


1 


15.7 


20.7 


15.6 


20.6 


14.5 


19.5 


2.22x10-1* 


2.53x10-16 


2 


22.3 


27.3 


22.0 


27.0 


20.6 


25.6 


1.34x10-18 


1.52x10-1* 


3 


28.9 


33.9 


28.3 


33.3 


26.7 


31.7 


8.21x10-18 


9.34x10-21 


4 


35.5 


40.4 


34.6 


49.6 


32.8 


37.8 


5.06x10-21 


5.76x10-23 



^Magnitudes are for bursts with all their fluence in Is. 
I'Values for HETE axe ergs/cm'^-s in the 5 — 7 eV band. 



Table 2. Instrument Summary 



Instrument Field of 


Exposure 


Sensitivity 


View (sr) 


Time 





Sonncberg'^ 


0.15 


40 m 


6 (ms)" 


Ondfejov'^ 


4.36 


4 hr 


3 (mv)'' 


Odessa'^ 


1.10 


30 m 


6 (mBT 


Dushanbe^ 


0.49 


1 hr 


6 (mv)^ 


ETC* 


0.75 


5 s 


11 (mv)^ 


GROCSE'' 


0.621 


0.4 s 


13 {mv)^ 


GTOTE^ 


0.06 - 0.12 


3.3-4 s 


12 - 14 (mv)^ 
8 X 10~^ ergs/cm^-s (5 


HETE^ 


1.7 


4 s 



'^Limiting magnitude for Is flash. 
''Sensitivity to field stars. 
^Greiner et al. (1994) 
'^Vanderspek et al. (1994) 

"Akerlof et al. (1994), values for the second generation camera are used 
here 

^S. Barthclmy (1995), private commimication (preliminary) 
sRicker et al. (1992); R. Vanderspek (1995), private communication 
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Table 3. Expected Brightness of Bursts for Exposure Duration 



Index 


ETC, GTOTE & GROCSE 


Ondfejov^ 


HETEP 


min(my) 


median(mv') 


min(mv') 


niedian(my) 


max{flux) 


median (/iux) 


a + 1 


9.9 


19.2 


18.2 


25.7 


1.57x10-12 


5.12x10-16 


-4 


-15.2 


-8.8 


-6.9 


-1.9 


4.75x10-4 


1.52x10-6 


-3 


-8.6 


-2.2 


-0.3 


4.7 


2.75x10"^^ 


8.82x10-3 


-2 


-2.0 


4.4 


6.3 


11.3 


1.60x10"^ 


5.14x10-11 


-1 


4.6 


11.0 


12.9 


17.9 


9.44x10-" 


3.03x10-13 





11.2 


17.6 


19.5 


24.5 


5.61x10"" 


1.80x10-15 


1/3 


13.4 


19.8 


18.7 


23.7 


1.02x10-13 


3.27x10-16 


1 


17.8 


24.2 


26.1 


31.1 


3.37x10-1^ 


1.08x10-1^ 


2 


24.4 


30.8 


32.7 


37.7 


2.04x10-1^ 


6.54x10-20 


3 


30.9 


37.4 


39.3 


44.3 


1.25x10-13 


4.00x10-22 


4 


37.5 


44.0 


45.9 


50.9 


7.68x10-22 


2.46x10-24 



lvalues for Sonneberg, Odessa, and Dushanbe are typically 2 magnitudes brighter than for Ondfejov. 
2 Values for HETE are in terms of ergs/cm2-s in the 5 — 7 eV band. 
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